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Using the optical system equation, the impact of a change in wavelength is 
determined. The rms phase variation at output aperture is critical and must be 
reduced by an amount equivalent to reduction in wavelength. If turret size is 
scaled with wavelength, the distortion due to the inviscid flow remains constant. 
Decreased wavelength offers the flexibility of greatly reduced turret weight, 
volume, and aerodynamic drag. Multiple turrets become feasible on an aircraft. 
The distortion due to the turbulent boundary layer becomes magnified by reduction 
of wavelength; however, if the turbulent boundary layer originates on the 
turret and if the turret is scaled, the boundary layer thickness decreases. 
Distortion becomes less significant in this case. 



iii 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

II. FAR-FIELD INTENSITY 1 

III. SCALING TURRET WEIGHT, VOLUME, AND DRAG 5 

IV. SCALING OPTICAL DISTORTION DUE TO INVISCID FLOW 6 

V. VISCOUS FLOW PHENOMENA LEADING TO OPTICAL DISTORTION 9 

A. Types of Viscous Flow Phenomena 9 

B. Turbulent Boundary Layers 10 

C. Scaling Propagation Through a Turbulent Boundary Layer of Fixed 

Thickness 14 

D. Scaling Propagation Through a Scaled Turbulent Boundary Layer . 15 

VI. SUI^IARY, CONCLUSIONS, AND RECOMflENDATIONS 17 

APPENDIX A 19 

REFERENCES 21 



ACKNOWLEDGJIENT 

The author wishes to thank his wife, Emily, for typing this report. 



iv 



LIST OF FIGURE AND TABLES 



FIGURE 

1. Diffraction angle for which the intensity is one half 



the maximum value 15 

TABLES 

I. Parameters for Turret Weight 5 

II. Change of Half Angle with Change in Wavelength 16 



V 



t 



AERO-OPTICS AT SHORTER WAVELENGTHS 



by 

Distinguished Professor Allen E. Fuhs 
Department of Aeronautics 
Department of Physics and Chemistry 
Naval Postgraduate School 
Monterey, California, 93940 



I . INTRODUCTION 

An extensive technological base for aero-optics has been developed by the 
laser community. The research and development work was oriented to the CO^ 
laser with wavelength of 10.6 microns. In recent times, lasers operating at 
much shorter wavelengths, e.g., iodine laser at 1.315 microns [1] or the 
excimer lasers in the blue end of the visible spectrum [2], have shown potential 
for superior performance. The information concerning laser turrets on aircraft 
needs to be reexamined for shorter wavelengths. 

Section II discusses the far-field intensity using an equation which was 
developed by Holmes and Avizonis [3]. The equation highlights the impact of a 
change in wavelength. Section III discusses scaling of turret weight, volume, and 
drag. Section IV discusses the scaling of phase distortion due to inviscid 
flow. Section V considers phase distortion to aircraft generated turbulence. 
Finally Section VI summarizes the paper and states conclusions relative to 
unknowns or areas of uncertainty. 

II. FAR-FIELD INTENSITY 

An equation for far-field intensity has been developed by Holmes and 
Avizonis [3]. The equation, which is useful for systems analysis, is 
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One additional term, D/FD^, has been added by Fuhs [4] to account for turbulence. 
The meaning of each term will be discussed now. 

The first term is the far-field intensity on axis for an Airy diffraction 
pattern. The mirror diameter is D; the laser power at the output aperture, 
i.e., laser telescope primary mirror, is P; the wavelength of the laser is A; 
and the range to the focal plane is Z. The shorter wavelength increases 
however, one may want to decrease mirror size, thereby decreasing turret weight 
and volume. Further, a decrease in mirror diameter, D, will decrease aerodynamic 
drag. Consequently, the ratio of D to A will be assumed to be constant, i.e., 

= A/D = constant (2) 



is a scaling factor. The turret will have a characteristic dimension L. An 
additional assumption is 

= L/D = constant (3) 



where is another scaling parameter. 
The term for jitter is 



1 . 
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(4) 



When is constant, the jitter term depends on the line-of -sight rms angular 
beam jitter, O. Many different phenomena influence the magnitude of O including 
the following: 

aerodynamic inputs to airframe 
aerodynamic input directly to the laser turret 
aerodynamic input due to control surface deflections 
unsteady aerodynamics in vicinity of laser turret 
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unsteady aerodynamics due to open port 
vibration due to aircraft components 
electrical noise within servo system 
laser turret mechanical jitter 

interaction of turret aerodynamics with remainder of aircraft 
Each item will not be discussed. When the laser turret shrinks in size to main- 
tain constant, the value of a must be maintained constant to avoid increased 
detrimental effects. 

From the causes of jitter listed above, the problem due to open port may be 

eliminated or decreased. With another wavelength, i.e., 1.3 microns instead of 

10.6 microns, and a much smaller size, a window may become feasible. 

The quantity is the loss due to the reflectivity of m-mirrors in the 

optical train. is the transmission of the window mentioned in the preceding 

paragraph. The quantity T^ is diffractive spillage around mirrors in the 

optical train. T, is a function of Fresnel number 
d 



D 
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where r is the radius of mirrors in the optical train, and £ is the spacing 
m 

between two mirrors. If one scales the primary mirror by ratio S^, then the 

diameter of a mirror in the optical train, D , undoubtedly will be scaled to 

m 

a smaller size. The consequence is increased spillage, i.e., a smaller value 

for T,. One assumes that £ remains fixed, 
d 

The transmission of the atmosphere is T and is given by 

a 



T = e 



( 6 ) 
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where y is an extinction coefficient. The extinction coefficient is the sum 
of four terms 

Y=a +a +k +k (7) 

' a s a s 

where is molecular absorption, is molecular scattering, k^ is aerosol 
absorption, and k^ is aerosol scattering. Each of the terms in equation (7) is 
dependent on wavelength. 

2 

The factor K is a consequence of nonuniform intensity, I, watts/cm , at 

the primary mirror of the laser turret. If is fixed, the intensity in the 

far-field does not change; see Appendix A. 

Aero-optics plays an important role in the factor 6 /A. Of course, 6 

^ ^ rms rms 

is the result of all optical aberrations starting within the cavity of the laser 
device. The aerodynamic window, the mirrors in the optical train, and the 
primary telescope mirror affect The beam propagates through a turbulent 

boundary layer and the variable density field generated by the inviscid flow over 
the turret. Both external flows contribute to 6 

rms 

Obviously, when A is decreased, 6 must be decreased by a like amount. 

rms 

Otherwise, is reduced. When A is reduced from 10.6 microns to 1.3 microns, 

i.e., by a factor of 8, 6 must be reduced also by a factor of 8. The most 

rms 

difficult aspect of scaling to shorter wavelengths will be to reduce 6 

rms 

A connection exists between in equation (A- 1) and 6 /A in equation (1) . 

rms 

For more details concerning the relationship between lb and 6 , see Born and 

rms 

Wolfe [5]. 

A term, D/FD^, appears in equation (1) to account for turbulence. The term is 
adequate for systems analysis but is not appropriate for a detailed examination of 
aero-optics. The question of turbulent flow is deferred to Section IV. 
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III. SCALING TURRET WEIGHT, VOLUME, AND DRAG 



Turret weight, W, can be scaled using the equation 

W2 D2 " 

— = r— 1 

Wi Id^J 



( 8 ) 



where n is a number depending on the turret optics and D is output aperture. 
Table I gives the value of n, W^, and for three turrets. 



Table I. Parameters for Turret Weight 



Type of Turret 


n 


, pounds 


, cm 


on-gimbal Cassegrain 


2.50 


1637 


60 


heliostat 


1.93 


1788 


60 


coelostat 


2.40 


2161 


60 



Changing from a CO^ laser to an iodine laser gives a wavelength ratio of 1/8. 
Assuming is fixed, the value of D can be reduced by a factor of 1/8; mirror 
diameter can be reduced from 60 cm to 7.5 cm. Use the numbers in equation (8) for 
an on-gimbal Cassegrain telescope. The result is 



W^ = 1637 




9.04 



The range of validity of the scaling law of equation (8) may have been exceeded; 
even so, the qualitative conclusion is that turret weight is very significantly 
reduced. 

Turret volume, V, follows a scaling law of 
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( 9 ) 



Once again, significant reductions occur. 
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Aerodynamic forces are given by an equation of the form 




( 10 ) 



where is a force coefficient. The drag is 




( 11 ) 



where is the drag coefficient. In equation (11), A is the cross sectional 
area of the turret exposed to the flow. Hence the drag varies as mirror diameter 
squared . 

When laser turrets weigh in tons, only one turret can be used on an aircraft. 
However, when the turret is reduced in weight and volume by a factor of 100 
or so, many turrets can be used. Instead of attempting to cover 2tt-s teradians 
with one turret, each of the small turrets can cover a much smaller solid angle. 
Furthermore, the turrets can be located at many different places on the aircraft. 
In the enthusiasm to add turrets, consideration must be given to fire control, 
pointing, and tracking. 

IV. SCALING OPTICAL DISTORTION DUE TO INVISCID FLOW 

Three lengths are of relevance to the scaling of optical distortion due to 
inviscid flow; these are beam diameter, D; wavelength, A; and characteristic 
turret dimension, L. One length is selected as a reference length to form 
ratios. Equations (2) and (3) indicate D is used as the reference. 

Compressible flow over a laser turret generates a spatially varying density 
field. A relation exists between refractive index and density 



n = 1 + (Kp/p^) 



( 12 ) 



The fractional changes are 



An _ ^ 



n p 



( 13 ) 
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By solving the equations of motion for flow over a turret, the density field can 
be obtained. Knowing the density field, the refractive index can be calculated. 

The optical path length, OPL, is defined as 

OPL = J n(x,y,z)ds (14) 

where s is distance along a ray. The optical path length is a geometrical 
length weighted by the local refractive index. The change in OPL, known as 
optical path difference, OPD, is 

6 = OPD = J n^(x,y,z)ds^ - y'n^(x,y,z)ds^ (15) 

where subscripts i and r indicate the i-th ray and a reference ray. The quantity 6 
is equal to the optical path difference. 

Define nondimens ional variables 

X = x’D 

y = y’D 

z = z ’ D 
s = s *D 

Substitution of the nondimensional variables into equation (14) gives 



OPL ^ ^ J" n(x* ,y’ ,z’)ds’ 



Likewise the optical path difference becomes 



6=D/n.dsl-D /nds’ 
J 1 1 J r V 



(16) 



(17) 



Hence OPL and 6 scale as mirror diameter D. One must show that the integrals 
in equations (16) and (17) are constant when turret size varies so that and 
are constant. 
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For blunt turrets, the equation for the velocity potential in compressible 



flow is 



1 3(1) 3^ 3^(1) ^ 3^d) 

2 

a 9x. 3x. 3x.9x, 9x.9x. 

1 k 1 k j j 



Introduce nondimensional variables 



D 



(18) 



a = a'a 

OO 



into equation (14). The result is 

^ 3 ^ ^ 3 ^({)* ^ 3 ^ 4 )' 

a'^ 3x1 3x,' 3xI3x,' 3x13x1 

1 K 1 Ic J J 



(19) 



A length scale parameter does not appear in equation (19) . Although not demonstrated 
here, the boundary conditions can be nondimens ionalized without a residual length 
scale. Consequently, the solution ({)* is independent of turret scale. 

The local Mach number is related to freestream Mach number by 



M = 



u u 



a'a 



= M M* 



( 20 ) 



The local Mach number determines the local refractive index. Since M* is 
independent of turret scale, the identical value of refractive index exists at 
corresponding points in two flows with different scale. Consequently, the integrals 
in equations (16) and (17) are identical for the two flows at different scale. 

The phase distortion, also known as optical path difference, relative to 
wavelength is 



6 

A 






D 

A 



/ 



n ds ' 
r r 



( 21 ) 
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Since the integrals are independent of scale, and since is fixed, the phase 

distortion does not vary. Likewise 6 /X does not vary with turret scale 

rms 

provided M , S, , and S are fixed. 

^ X z 

For small perturbation turrets, a similar argument can be made. 

V. VISCOUS FLOW PHENOMENA LEADING TO OPTICAL DISTORTION 

A. Types of Viscous Flow Phenomena 

Viscosity dominates and determines, at least in part, the gross 
characteristics of several different types of flow. These include boundary 
layers, shear layers, open port cavities, wakes, and shock waves. Although not 
a flow of the same nature as those listed above, jet engine exhausts should be 
included; exhausts can have a detrimental effect on beam propagation. At the 
speeds and geometric scales involved with aircraft, the viscous flows lead to 
turbulence. Extensive literature [6-21] exists on the propagation of electro- 
magnetic waves through a turbulent medium. 

Boundary layers surround the aircraft surfaces. The boundary layer on 
leading edges and forward surfaces is laminar with transition occurring at a 
critical Reynolds number. Flow separation leads to shear layers. Shear layers 
separate two flow regions with different flow velocities. 

Due to the extremely high irradiance in the laser beam, the turret may have an 
open port. Windows are not available which can survive the high irradiance. The 
interior of the turret is a cavity in which unsteady, turbulent flow exists. The 
flow in a cavity may have regularity in the sense that discrete vortices can be 
identified . 

Wakes are turbulent and adversely affect propagation. When selecting a 
location for a laser turret, any location within or near a wake should be avoided. 
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Shock waves create large changes in density and refractive index. Shock 
waves which originate or reflect from boundary layers create intense turbulence. 
Depending on conditions, a triple point may occur in reflection of an oblique 
shock. Downstream of the triple point, a shear layer exists. Also downstream of 
a reflected shock, a separated flow region may occur. 

Two cases can be identified. In one case, the laser beam propagates through 
a region disturbed by a shock wave. In the other case, the shock wave originates 
or reflects near the aperture of the turret. The latter case causes the greater 
difficulty with optical distortion. 

B. Turbulent Boundary Layers 

The boundary layer may originate on aircraft surfaces far upstream of the 
aperture. If this is the case, then the boundary layer will be thick. When 
the laser turret protrudes into the flow, the boundary layer may begin on the 
turret itself. In that case, the boundary layer is thin. One can state 
unequivocally that the thinner the boundary layer, the lesser the adverse 
influence on beam quality. 

In addition to the distance which the boundary layer has traversed, other 
factors influence the boundary layer. Mach number characterizes the compressibility 
effects. At higher Mach numbers, the boundary layer tends to grow in thickness. 

Wall curvature, which usually is accompanied by pressure gradients, has an 
influence on boundary layer thickness and velocity profile. When the wall is 
concave to the external flov;, Gortler vortices may occur. Cooling or heating can 
alter boundary layer properties. In fact, cooling may be one approach to flow 
control. Surface roughness may amplify turbulent intensities. Surface roughness 
may be due to rivet heads, structural joints, insects which have smeared the sur- 
face, etc. 
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Pressure gradients have been mentioned earlier. A gradient in which the 
pressure increases in the flow direction is termed an adverse gradient. 

Conversely, a gradient with decreasing pressure in the flow direction is a 
favorable gradient. Adverse pressure gradients lead to flow separation as well 
as thicker boundary layers. Favorable pressure gradients delay transition from 
laminar to turbulent flow. Favorable pressure gradients tend to cause thinner 
boundary layers. A boundary layer in a favorable pressure gradient rarely separates. 

Shock wave, boundary layer interaction has been mentioned. If the inter- 
action occurs at the aperture or immediately upstream of the aperture, severe beam 
degradation is likely. 

For turbulent incompressible boundary layers, the thickness, 6, is given 
by [22] 

- = 0.37 Re (22) 

X X 

According to equation (22), 6 increases as x to *f 0.8 power. For a laminar, 
compressible boundary layer, the thickness is given by [23] 



6 = [1 + (23) 

where and are constants. Turbulent compressible boundary layers show 
similar trends with increasing Mach number. 

A boundary layer in compressible flow has a temperature gradient across 
the boundary layer. The gradient in temperature is due to viscous dissipation. 

A quantity known as the recovery factor r is defined for an adiabatic wall; r is 



r 




Pr 



(24) 
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The temperature of the wall is T^, the temperature of the air at the outer edge 
of the boundary layer is T^, and the stagnation temperature is . Equation (24) 
shows r equal to the Prandtl number, Pr. An effective Prandtl number can be 
defined for a turbulent boundary layer and will have a value near unity. 

Hence the wall temperature will be equal to the stagnation temperature when 
Pr = 1.0. 

The rms density fluctuations in the turbulent boundary layer are related to 

the difference in density at the wall, p , and at the edge of the boundary 

w 

layer, p^. According to Sutton [24] 



p = 0.1(p - p ) 

rms oo w 



(25) 



for boundary layers with supersonic external flow. For subsonic boundary layers, 
the constant is larger than 0.1 shown in equation (25). The density at the wall 
is given by 



P,, = 



w RT 



(26) 



w 



The rms fluctuation in refractive index is related to density by 



rms rms 

= K 

n p 



(27) 



where K is the constant from equation (12). 

An attenuation coefficient, a, can be defined for propagation through a 
layer of a turbulent medium. The intensity ratio is 



I -a6 



(28) 



where the layer thickness is 6, and is the intensity of the beam entering the 
layer. The attenuation is due to scattering by the turbulent boundary layer. 
Sutton [25] gives a formula for a as follows: 



12 



2 2 

a = 2k n A 
rms 

where k is the wavenumber equal to 2 tt/X and A 
The macroscale is obtained from 



( 29 ) 

is the macroscale of the turbulence. 



A = 



y' C(r)dr 
0 



(30) 



where C(r) is the correlation function for refractive index in an isotropic 
turbulent medium. The correlation function is given by 

cx> 

J An(r + R)An(r)47rR^dR 

C(r) = (31) 

J [An(R) l^AirR^dR 
0 

where An is the deviation from the average, n, 

An = n(r) - n (32) 

Sutton [25] reports values of the scale of turbulence equal to about 0,1 6, i,e,, 
the scale of turbulence is approximately 1/10 of the boundary layer thickness. 

For propagation through turbulent boundary layers, several characteristic 
lengths occur as follows: 

A macroscale of turbulence 

6 boundary layer thickness 

D aperture for laser beam 

A wavelength of laser radiation 

Many regions for turbulent propagation can be characterized by the length scales 
along with the extinction number a6. 
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C. Scaling Propagation Through a Turbulent Boundary Layer of Fixed Thickness 

Two parameters can be used to characterize the degradation in beam quality 

2 

due to turbulence. One is the decrease of intensity, watts/cm , on the axis 
of the beam in the far-field. For large values of D/A, equation (28) provides a 
means to calculate I/I^. The second parameter is the half-intensity angle for 
the beam. The half-intensity angle, defines a cone with a vertex angle 

equal to surface of the cone, the intensity is one-half of the peak 

intensity. 

Figure 1 shows the half-intensity angle as a function of the ratio D/A. The 
ordinate is the variable . This variable compares the half- 

intensity angle for the scattered beam with the half-intensity angle for a 
diffraction limited beam. 

Consider a CO 2 laser beam at 10.6 microns from a D = 60 cm mirror. Assume the 
boundary layer thickness, 6, is 10 cm, and the macroscale. A, is 1 cm. For 
M = 1.5 flow with turbulent Prandtl number of unity and p^ = 1 atmosphere, 

= 1.23 kg/m^, and T = 290°K, the values of p and p are 0.848 kg/m^ and 

^00 » 00 » '^rms ^ 

3 -4 -6 

0.038 kg/ra . Using K = 2.4 x 10 , the value of n /n is 7.44 x 10 . The 

rms 

extinction coefficient, which is calculated using equation (29), is 0.39 per meter. 
The extinction number, a6, is 0.039, and D/A is 60. The corresponding point can 
be plotted in Figure 1, giving almost on the dotted line. 

Consider now an iodine laser at 1.315 microns with a turret for which 
D/X = has been maintained constant. Hence D = 7.44 cm; also, D/A is 7.44. 

The extinction coefficient grows to 

a = (0.39)(8.06)^ = 25.34 
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Figure 1. Diffraction angle for which the intensity 
is one half the maximum value. Reproduced 
from Sutton [ 25 ] . 



and the extinction number is 2.53. The point for the iodine laser can be plotted 
in Figure 1, giving a half angle parameter of about 0.6. The increase in half 
angle is about 20 per cent. 

For a thicker turbulent boundary layer with 6 = 20 cm, the results are more 
dramatic as shown in Table II. The half angle parameter increases from 0.48 to 
1.2, which is a factor of 2.5. 

It should be emphasized that equation (28) does not apply when D/A is less 
than approximately 10. The attenuation is less than that predicted by equation (28) 
when D/A < 10. 

D. Scaling Propagation Through a Scaled Turbulent Boundary Layer 

This subsection discusses the case for which the turbulent boundary layer 
originates on the laser turret instead of an aircraft surface. From equation (22), 
the ratio of boundary layer thicknesses is 
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Table II. Change of Half Angle with Change in Wavelength 




wavelength, A, micron 10,6 


i.3i; 


boundary layer thickness, 6, era 20 


20 


macroscale. A, era 2 


2 


extinction coefficient, a, per meter 0.78 


50.6 7 


extinction number, a6 0.156 


10.13 


ratio D/A 30 


3.72 


ratio D/A 56,603 


56,603 


Mach number M 1.5 


1.5 


pressure, p^, atmosphere 1 


1 


3 

density, p^, kg/m 1.23 


1.23 


temperature, T^, 290 


290 



0.8 , 0.8 

6 X L 

^2 ""l h 

Further, the ratio of macroscales of turbulence is 


(33) 


Ai - 

The extinction number scales as 


(34) 


aS) A A L L L., 

Equation (35) assumes is not changed. The ratio D/A scales as 


(35) 
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0.8 



0.2 



D/A) 2 

D/A)j^ 



= ^ ^ A. 

- Li A2 - 4 ^L2^ 



Using equations (35) and (36) along with data of 
are found for A = 1,315 microns: 




Table II, 



( 36 ) 

the following values 



D/A = 19.8 
a6 = 1.89 



Plotting the preceding points in Figure 1 yields a ^bout 8 per cent larger. 

The impact of reduced A is much less when the boundary layer originates on the 
turret. 



VI, SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 

The analysis has assumed that S^^ = A/D and S^ = L/D were maintained constant 
as A was decreased. The conclusions are valid for this case. 

Using the optical system equation of Holmes and Avizonis [3], various conse- 
quences of reducing A were examined. In summary, the following should be noted: 

6 Phase distortion must be reduced by a factor of 8 for a change 

rms ^ ^ 

from CO^ to iodine laser. 

K A factor which depends on the spatial distribution of intensity 

at the output aperture will not vary, 

T^ Diffractive spillage will tend to increase reducing the value of 

O rms jitter can be maintained constant without change in l^^* 

Detailed comments were not made about atmospheric extinction coefficient; 
however, scattering tends to increase as A decreases. 

Turret weight and volume as well as turret aerodynamic drag will be reduced 
dramatically. Much greater flexibility in the number and location of turrets is 
a consequence. Further, the possibility of obtaining a suitable window should 
improve; open port problems are thereby circumvented. 
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The analysis shows that the phase distortion 6/A for inviscid flow is givei 



by 



6 

A 



(constant) 



D 

A 



( 37 ) 



Since D/A is constant, the phase distortion due to inviscid flow is independent 
of wavelength. The statement is true provided and are fixed. Nonetheless 

severe degradation due to the inviscid flow field can occur; hence, the problem 
cannot be ignored. 

For turbulent boundary layers, the value of the extinction number, a6, 
increases as A decreases. Also, the ratio D/A decreases as A decreases. Since 
extinction and scattering depend on a6 and D/A, the trends are not favorable. 

The half-intensity angle increases, and extinction increases as A decreases. 
Turbulent boundary layers require more attention. 

Any change which causes a thinner turbulent boundary layer tends to improve 
beam quality. Mounting the laser turret at points farther forward decreases 6. 
Cooling the wall in the supersonic case also should be an improvement- Finally, 
attention to surface roughness should help reduce scattering and attenuation. 

When the turbulent boundary layer originates on the laser turret, the 
boundary layer thickness decreases as the turret decreases in size. Hence, the 
impact of the turbulent boundary layer on almost unchanged by a change 

in A; however, serious degradation in intensity due to a6 remains. 
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APPENDIX A 



Consider the change in intensity at some point P in the far-field of an 
aperture illuminated in phase and amplitude by 

F(C,n) = (A-l) 

Silver [26] on page 173 shows the amplitude and phase of the wave at 
point P is given by 



= 



AR 



-jkR 



/ 



F(C»D)exp[jk sin 6(^ sin (}) + n cos (j)) ]d^dn 



(A-2) 



Intensity at point P is 



P P 



(A-3) 



Introduce nondimensional variables 

c = C*D n = n’D (A-4) 

where D is the diameter of the aperture. The integration in equation (A-2) 

2 

extends over A = ttD /4. Examine first the exponent in the integrand of 
equation (A-2) ; it becomes 

jk sin 0[C sin (}) + n cos (j)] = jkD sin 6[C* sin 4) + ri ’ cos (j)] (A-5) 

Note that 



kD = 2i\/S^ (A-6) 

The complex function F(C,n) is exactly the same as equation (A-l) except primes 
are added to ^ and n* The integral in equation (A-2) , which is rep>resented by 1^, 
becomes 



D^I^ = J* F(C’ ,n’)exp[ j (2Tr/S^)sin 0(5’ sin (}) + n’ cos (}))i]d5’dn* (A-7) 

A’ 

The intensity at point P is obtained by combining equations (A-2), (A-3) , and (A-7) 
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(A-8) 



2 2 2 
D I 

n _ n 

” \2 2 ~ 2 2 
S7R 



Consider two cases 1 and 2 which differ in the wavelength of radiation but for 



which is constant. The ratio of intensities becomes 

2 2 

!p2 

^Pi ■ r2 
nl 2 



(A-9) 



In equation (A-9) , the integrals, I , cancel when P is located at the same 



point specified by R,0,(}). Since R^ = R^, the ratio of intensities is identical. 
Consequently, the factor K in equation (1) is identical for cases 1 and 2. 
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